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High resolution cosmological N-body simulations of four galaxy-scale dark matter halos are com- 
pared to corresponding N-body /hydrodynamical simulations containing dark matter, stars and gas. 
The simulations without baryons share features with others described in the literature in that the 
dark matter density slope continuously decreases towards the center, with a density pdm oc r ~ 13±0 - 2 ] 
at about 1% of the virial radius for our Milky Way sized galaxies. The central cusps in the simula- 
tions which also contain baryons steepen significantly, to pdm oc r _1 - 9±0 ' 2 i with an indication of the 
inner logarithmic slope converging. Models of adiabatic contraction of dark matter halos due to the 
central build-up of stellar/gaseous galaxies are examined. The simplest and most commonly used 
model, by Blumenthal et at, is shown to overestimate the central dark matter density considerably. 
A modified model proposed by Gnedin et al. is tested and it is shown that while it is a considerable 
improvement it is not perfect. Moreover it is found that the contraction parameters in their model 
not only depend on the orbital structure of the dark-matter-only halos but also on the stellar feed- 
back prescription which is most relevant for the baryonic distribution. Implications for dark matter 
annihilation at the galactic center are discussed and it is found that although our simulations show 
a considerable reduced halo contraction as compared to the Blumenthal et al. model, the fluxes 
from dark matter annihilation is still expected to be enhanced by at least a factor of a hundred 
as compared to dark-matter-only halos. Finally, it is shown that while dark-matter-only halos are 
typically prolate, the dark matter halos containing baryons are mildly oblate with minor-to-major 
axis ratios of c/a = 0.73 ± 0.11, with their flattening aligned with the central baryonic disks. 

PACS numbers: 95.35.+d, 98.62.Gq, 95.10.Ce, 95.30.Lz 



I. INTRODUCTION 

There are still a multitude of open questions regarding 
the distribution of dark matter in galactic halos, the ef- 
fect of baryons upon the structure of dark matter halos 
being one. Many simulations of galactic halos only take 
into account the dark matter (e.g. [3, S, H, 0> S @| ) i but 
one would expect the baryonic component of the galaxy 
to behave very differently since it is able to cool (dissi- 
pate energy), and contract considerably. This is indeed 
observed in simulations containing baryons and also in 
nature, where the baryons form a disk and/or bulge at 
the center of apparently much more extended dark mat- 
ter halos. 

It has long been realized that this ability of baryons 
to sink to the center of galaxies would create an en- 
hanced gravitational potential well within which dark 
matter will congregate, increasing the dark matter den- 
sity there. To model this effect it is common to use 
adiabatic invariants or some small modification of them 
SlllUlliSll. Such models are fre- 
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quently used and it is of particular interest to test the 
validity of these models [IH, [l?], HI] currently because 
of recent advances in mapping the velocity field of the 
Milky Way [l|| and other galaxies. Furthermore, up- 
coming gamma ray experiments will look for the flux due 
to the self-annihilation of weakly interacting dark matter 
candidates from, e.g., the galactic center. Since this flux 
is proportional to the dark matter density squared, pre- 
dictions for, and conclusions from, the data will depend 
strongly on the details of the effect of baryonic pinching 
upon dark matter halos. 

In this work we aim to investigate the effects of baryons 
on dark matter halos in disk galaxies and test the most 
common models of adiabatic contraction by comparing 
recent cosmological N-body/hydrodynamical simulations 
of galaxies containing dark matter and baryons to results 
from N-body dark-matter-only simulations of the same 
halos. Fully cosmological simulations starting at high 
redshifts are used. The simulation results are known from 
previous studies to produce overall realistic gas and star 
structures for spiral galaxies 0, [U H3| , even though the 
numerical resolution is still far from being able to resolve 
any small scale features observed in real galaxies. The 
most important dynamical property in this paper is the 
creation of stable disk and bulge structures both for the 
gas and star components. The angular momentum prob- 
lem is overcome by stellar feedback, implying that the 
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matter is not too centrally concentrated; a gene ric prob- 
lem of early galaxy simulations (see, e.g., [23| and ref- 
erences therein on forming disk galaxies in simulations). 
With the baryonic disks and bulges formed fully dynam- 
ically the surrounding dark matter halo response should 
also be realistically predicted inside the simulated spiral 
galaxies. 

The paper is organized as follows. In section |TT] we 
present the simulations, and then focus in section IIII1 
in particular, upon profile fits to the density of the 
dark matter halos in the simulations with and without 
baryons. In section ITVl we investigate different prescrip- 
tions which aim to predict the effects of baryons upon 
dark matter profiles by testing whether they are able 
to reproduce the dark matter profiles observed in the 
simulations with baryons. We then comment in section 
fVl on how these results might change the expected flux 
of gamma rays from dark matter annihilation in spiral 
galaxies, including the Milky Way. Finally, we analyze 
in section I VII the nonsphcricity of the different compo- 
nents of the simulations - dark matter, gas and stars - 
to attempt to further quantify their effects upon each 
other before we summarize our results in section IVIII 



II. SIMULATIONS 

The simulated galaxies used in our work consist of two 
Milky Way sized galaxies with virial radii (at z=Q) of 
r2oo ~ 200 kpc and two smaller galaxies with virial radii 
of around 100 kpc. The two larger galaxies are labeled 1 
and 2 and the two smaller galaxies 3 and 4; the number 
simply being an identifying label. Here we have followed 
common practice and defined r2oo as the radius of the 
sphere enclosing the mass M200 within which the mean 
density is 200 times the critical density, p c — 3Hq/8itG. 

All four galaxies have been extracted from fully cos- 
mological simulations using the Hydra code and an im- 
proved version of the Smoothed Particle Hydrodynam- 
ics code TreeSPH. Since the software generating the 
disk galaxies has been used in many previous works (see, 
e.g.,[2fj, |2l|, |22| ), we will only briefly mention the main 
features of the numerical code. 

The simulations are performed in a AC DAI cosmology 
with tt M = 0.3, fl A = 0.7, H = 100/ikm s^Mpc" 1 = 
65km s _1 Mpc _1 and with the matter power spectra nor- 
malized such that the present linear root mean square 
(rms) amplitude of mass fluctuations inside 8/i _1 Mpc is 
erg = 1.0. There is still some uncertainty in the accepted 
values of as and h, but we have tested that none of the 
general conclusions obtained in this paper are affected by 
small changes in these parameters. 

The galaxies are generated by first performing a dark- 
matter-only simulation, using the Hydra code, with 
128 3 particles in a box of comoving length of 10 hT 1 
Mpc and starting at redshift Zi = 39. After running this 
simulation, galactic size objects are identified. The simu- 
lations which include baryons (or alternatively with only 



dark matter) are then set up, particles within 4r v ; r at 
z = are traced back to their initial conditions, the dark 
matter particle mass resolution is increased by up to a 
factor of 64, and one SPH particle, i.e. baryonic mat- 
ter, per dark matter particle is added (keeping the total 
mass and fixing the baryonic fraction / to 0.15). The 
simulations are then rerun with the improved TreeSPH 
code; incorporating star formation, stellar feedback pro- 
cesses, radiative cooling and heating, etc. The final result 
are qualitatively similar to observed disk and elliptical 
galaxies at z = 0, a result which is mainly possible by 
overcoming the angular momentum problem by an early 
epoch of strong, stellar energy feedback in form of SNII 
energy being fed back to the intrastellar medium. 

To study the effect of baryons on dark matter ha- 
los in Milky Way like galaxies we use four simulated 
disk galaxies with the highest available resolution at 
our disposal. By comparing simulations with different 
resolutions we infer that the results are robust down 
to an inner radius r m j n , which is 2 times the gravita- 
tional softening length. This is in approximate agree- 
ment also with other commonly used convergence cri- 
teria, such as r m in ps A^c/^oo (see, e.g., [1, HJ] and 
references therein). For the larger galaxies 1 & 2 we de- 
duce r m i n = 2kpc whereas for the two smaller 3 & 4 
^min = lkpc. The spiral galaxies containing dark matter 
and baryons will be labeled with 'S' (e.g. SI) whereas 
those containing only dark matter will be labeled with 
'DM' (e.g. DM3). There are some notable differences be- 
tween the four galaxies, for example the large galaxy SI 
has a very pronounced flat gas and stellar disk with a star 
bulge, whereas galaxy S2 is strongly barred. The gas in 
the two smaller galaxies has a very definite disk structure 
and the stars exhibit both disks and central bulges which 
are more centrally concentrated than the ones in the two 
larger galaxies. 

A summary of the main parameters of the simulated 
galaxies is given in Table [J (see [20L l2ll . [22j for further 
details). 



III. DARK MATTER HALO RADIAL PROFILES 

In the section on nonsphcricity we will address the non- 
spherical aspects of the simulations in more detail, in 
particular the effect of the triaxiality of the dark mat- 
ter halo upon the baryons and vice-versa. However be- 
fore that, in order to directly compare the profiles with 
most others in the literature, we assume spherical sym- 
metry and fit the galaxies with common radial profiles. 
One parametrization which assumes two asymptotic ra- 
dial power law behaviors at both small (7) and large (/?) 
radii is known as the 'a/37' profile (or the Zhao profile), 
where the density as a function of radius is given by the 
expression 



(r/r.)T [1 + (r/r.)«] ^ 
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TABLE I: Main properties 


of dark matter halos 


and g 


llaxies at 


z=0 








Simulation 


SI 


DM1 


C O 

bl 


DMz 


C Q 

BO 


Divio 


o4 


DM4 


Virial radius r2oo [kpc] 


209 


211 


200 


201 


100 


102 


98.3 


97.5 


Total mass M 20 o [10 n M Q ] 


8.9 


9.3 


7.8 


8.0 


1.0 


1.1 


0.93 


0.91 


Number of particles -/V200 [xlO 5 ] 


3.6 


1.2 


3.5 


1.0 


3.2 


0.98 


3.1 


0.90 


DM particle mass mj m [10 6 Mq] 


6.5 


7.6 


6.5 


7.6 


0.81 


0.95 


0.81 


0.95 


SPH particle mass mbaryon [1O 6 M0] 


1.1 




1.1 




0.14 




0.14 




Grav. soft length DM ed m [kpc] 


1.0 


1.0 


1.0 


1.0 


0.5 


0.5 


0.5 


0.5 


Grav. soft length SPH ebaryon [kpc] 


0.6 




0.6 




0.3 




0.3 




Characteristic circular speed V c [km/s] a 


245 




233 




124 




122 




Specific angular momentum (bulge+disk), j„ [kpc km/s 


447 




303 




144 




153 




Specific angular momentum (cold gas), j cg [kpc km/s] 


1895 




1055 




1005 




1093 




Star formation rate (SFR) [M@ /yr] 


1.4 




1.7 




0.13 




0.15 




b = SFR/(SFR) 


0.23 




0.32 




0.13 




0.16 




Baryonic disk + bulge mass [1O 1O M0] 


7.17 




5.79 




1.13 




0.98 




Baryonic Bulge-to-disk mass ratio 


0.19 




0.80 




0.76 




0.60 





a V c is determined as in 120 



where a governs the radial rate at which the profile in- 
terpolates between the asymptotic powers —7 and —(3. 

Often in the literature various constraints on these pa- 
rameters are assumed, resulting in subclasses of less gen- 
eral profiles. In particular the density profile in the outer 
part of the galaxy halo, f3, is often assumed to be 3, 
which is the canonical value found for halos in simula- 
tions invoking only dark matter. For this reason many 
authors use (a, /3,7) = (1,3,7) an d is commonly known 
as the generalized Navarro, Frenk and White (NFW) pro- 
file (with the standard NFW profile having 7 = 1) 
We in general find a significantly better fit to our profiles 
if we also leave a and j3 as free parameters. Moreover, it 
is obviously not certain that profiles used to character- 
ize the halos of dark-matter-only simulations should also 
give good fits to the dark matter halos formed in simu- 
lations also containing baryons. We have therefore done 
least x 2 fits to the profiles leaving the four parameters 
in Eq. (UJ free (po we always constrain by the total mass 
in the fitting range which we set to be between r mm and 
?"20o)- The profile fits can be seen in Fig. [TJ [2] as solid 
lines and the parameter values are given in Table HT1 

It is useful to keep in mind that with four free param- 
eters in the 'a/97' profile, there are some degeneracies in 
the inferred parameter values (25j . The numbers given 
in Table []TJ aim to give a good parametrization of the 
density profiles in the fitted range, and do not neccesar- 
ily claim to represent profiles that could be extrapolated 
into smaller radii with confidence. 

We have also fitted the dark matter halo density distri- 
butions over the same radial range using exponential pro- 
files, where the logarithmic slope changes continuously 
with radius as recently suggested by Navarro, Frenk and 



White 



p(r) 



P-2 exp 



r 
r-2 



1 



(2) 



In this profile, p_2 and r_2 correspond to the density and 
radius where p oc r~ 2 . The best fit values can be found 
in Table IIIII and the profile fits can be seen in Fig. [TJ [2] 
as dashed lines. 

Finally, we have fitted single power laws to the cen- 
tral dark matter density profile (r min < r < 0.05r 2 oo) to 
determine the averaged logarithmic slope of the resolved 
central cusp. These slope values are found in Table [TV] 
This can also partly be compared to the central asymp- 
totic logarithmic slope in Eq. (JJ, which is —7 (whereas 
in the profile in ^ the logarithmic slope is continuously 
decreasing towards the center). 

To describe the dark matter halos we actually fit 
dm/dr = Airpr 2 since this is representative of the actual 
data in the simulation. Logarithmic binning of the ra- 
dius is used to minimize the effect of substructure, which 
becomes more important at large radii, while at the same 
time capturing the behavior of density profile in the cen- 
tral part. 

In order to obtain an estimate of the variance with 
which to do the fits and for our subsequent analyses 
throughout the paper, we take five snapshots of each sim- 
ulation at different times. These snapshots correspond to 
today's epoch and to four successive earlier times with 
At = 200 Myr, and from these we find their one stan- 
dard deviation (unbiased) mean square dispersion. This 
At is long enough so that the particles in the inner re- 
gions of the simulations will have had time to completely 
change their positions with respect to the center of the 
galaxy, but a short enough timescale in total so that the 
overall mass profile of the halo has changed very little 
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FIG. 1: Top panels: dark matter density profiles from galaxy simulations including baryons; galaxy SI (left) and S2 (right). 
Bottom panels: density profiles from halo simulations including only dark matter; halo DM1 (left) and DM2 (right). The 
solid, dashed and dot-dashed curves show the best fit from the parametrization given in Eq. fTp, f^) and a single power law, 
respectively. Long arrows show the lower resolution limit (r m i n ) and the virial radius (V200 ), respectively. The shorter arrows 
indicate the upper limit for the single power law fits ( 0. 05 r2oo ) ■ 



TABLE II: Best fit parameters for model Eq. £Jp to the spherical symmetrized dark matter halos. 



Galaxy sim. with(without) baryons 


r s [kpc] 


a 





7 


Xdof [46 dof] 


SI (DM1) 


44.9 (36.0) 


1.76 (1.47) 


3.31 (3.36) 


1.83 (1.42) 


1.4 (1.0) 


S2 (DM2) 


150 (85.1) 


0.486 (0.588) 


4.21 (4.79) 


1.49 (0.850) 


1.5 (1.2) 


S3 (DM3) 


13.3 (14.3) 


12.4 (1.387) 


2.48 (2.74) 


2.07 (1.70) 


1.0 (2.1) 


S4 (DM4) 


56.1 (10.3) 


2.75 (0.915) 


3.48 (3.00) 


2.20 (1.36) 


2.1 (1.0) 



(for all galaxies, the dark matter mass within 50 kpc in- 
creases by less than about 1% over the total period of 1 
Gyr). 

This method of estimating the variance in each bin has 
the effect of suppressing the influence from temporary 
small scale inhomogeneities in density due to dark mat- 



ter subhalos (which are most common at large radii) and 
at the same time retaining the Poisson population vari- 
ance in other bins. It turns out that the variances are 
approximately equal to the Poissonian values for most 
radial bins. This will serve as the estimate of the uncer- 
tainties in our subsequent investigations. A more detailed 
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FIG. 2: Top panels: dark matter density profiles from galaxy simulations including baryons; halo S3 (left) and S4 (right). 
Bottom panels: density profiles from halo simulations including only dark matter; halo DM3 (left) and DM4 (right). The 
solid, dashed and dot-dashed curves show the best fit from the parametrization given in Eq. {7p, |]|) and a single power law, 
respectively. Long arrows shows the lower resolution limit (V m i n ,) and the virial radius (V200 ), respectively. The shorter arrows 
indicate the upper limit for the single power law fits ( 0. 05 r2oo ) ■ 



TABLE III: Best fit parameters for model Eq. £§|) to the spherical symmetrized dark matter halos. 



Galaxy sim. with(without) baryons 


r-2[kpc] 


a 


Xdof [48 dof] 


SI (DM1) 


11.9 (18.5) 


0.185 (0.247) 


3.1 (1.5) 


S2 (DM2) 


7.66 (18.7) 


0.117 (0.226) 


1.5 (1.5) 


S3 (DM3) 


2.30 (6.70) 


0.0728 (0.132) 


2.6 (2.1) 


S4 (DM4) 


0.750 (6.31) 


0.0507 (0.153) 


3.2 (0.91) 



analysis of the actual uncertainties which would take into 
account, e.g. systematic radial dependencies from numer- 
ical and resolution effects or even effects from the imple- 
mentation of physical processes themselves, is very dif- 
ficult to achieve and is beyond the scope of this paper. 
From the five time frames no strong correlations between 
our bins were found and we will in practice not take into 



account such eventual correlations. 

Tables ILTlirVI contain significant information: For the 
'a/37' profiles one has p oc r^ 13 at large radii, and even 
when (3 is left as a free parameter, values for [3 around 3 
emerge, albeit with a rather large scatter. Moreover, the 
X 2 per degree of freedom is significantly smaller when us- 
ing the 'a/?7' profiles rather than the exponential profiles 
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FIG. 3: Logarithmic slope of the dark matter density profiles 
for all halos, plotted versus radius. Solid/red curves from the 
simulation including baryons, and dashed/blue curves from 
the simulations including only dark matter. Left panel: the 
larger, Milky Way sized, galaxies (1 & 2). Right panel: the 
two smaller galaxies (3 & 4)- 



TABLE IV: Best fit parameters for single power law p oc r 1 
fits to the central dark matter density profiles (r m in < r < 
0.05r 20 o ). 



Galaxy sim. 
with(without) baryons 


7 


Xdof 

[49 dof] 


SI (DM1) 


1.91 (1.56) 


0.65 (1.2) 


S2 (DM2) 


1.91 (1.41) 


0.76 (1.1) 


S3 (DM3) 


2.13 (1.79) 


0.90 (0.92) 


S4 (DM4) 


2.21 (1.68) 


1.6 (0.80) 



for the simulation with baryons. 

From the best fit parameters in Table [ill it follows 
that the asymptotic central logarithmic slope (-7) for the 
dark-matter-only simulations average to -1.1 for the two 
larger halos DM1 and DM2 and a somewhat steeper slope 
of -1.5 for the two smaller galaxies. These results are, to 
the inner resolved radii, in agreement with other recent 
simulations including only dark matter (see e.g. [26[ and 
references therein). 

The effect that the presence of baryons has upon the 
central slope of the dark matter density profile is very 
pronounced, pushing the average asymptotical central 
logarithmic slope up to -1.7 for the larger galaxies and to 
-2.1 for the two smaller galaxies. Another way of seeing 
this systematic steepening of the profiles due to the pres- 
ence of baryons is to look at the fits to the exponential 
profile ([2]) presented in Table IIIIl For all four galaxies, 
the presence of baryons brings in the radius at which the 
density is dropping as p oc r~ 2 (for the smaller galaxies, 
this radius becomes comparable to r m i n ). 

Since all simulations have limited resolution one should 



be cautious in extrapolating profiles inside the re- 
solved radius of the simulations. To explore the inner 
slope further we show in Fig. [3] the logarithmic slope 
dlogp/dlogr. Note from this figure that there is no 
sign of convergence of the central slope in the simula- 
tions without baryons, whereas in the simulation includ- 
ing baryons the slope change is drastically less and might 
already have converged to a defined value. To calculate 
the logarithmic slope profiles without excessive particle 
noise we average both over the five time frames and over 
the five nearest radial neighbor bins (corresponding to 
logarithmic smearing window of 20% of our logarithmic 
radius range r min < r < r^oo)- This eliminates most of 
the fluctuations without biasing the slope significantly. 

For the two larger galaxies Fig. [3] clearly demon- 
strates that the logarithmic slope is always continuously 
changing for the dark-matter-only simulations, while for 
the simulations including baryons the logarithmic slope 
derivative is drastically less in the central regions. This 
may indicate that the inner logarithmic slope converges 
to a value close to -1.9, i.e. close to an isothermal sphere. 
The picture is slightly less clear for the two smaller galax- 
ies, but there the logarithmic slope also flattens out in 
the inner region to a value roughly around -2.1 (although 
the bin-to-bin scatter is somewhat too large to draw any 
firm conclusions). 

If we only fit the central part of the dark matter density 
{^min < r < 0.05r2oo) a single power law fit should work 
well for the simulation including baryons and we deduce 
from Table ITVl a inner logarithmic slope close to -1.9 for 
the two larger galaxies and an average inner slope of -2.2 
for the two smaller galaxies. 

Galaxy S2 is different compared to the other galax- 
ies in several ways. For instance, for the dark-matter- 
only simulation of the halo (DM2), the best fit to the 
'a/37' density profile yields a asymptotic central loga- 
rithmic slope of only —7 = —0.85 and a very large scale 
radius r s . One thing that separate this galaxy from the 
others in its evolution is that it has experienced a late 
time merger and that there is not as strong disk struc- 
ture, but rather a pronounced bar structure in the stars. 
Despite this the dark matter density profile is not entirely 
different, which seems to support the results reported in 
[27J - where they find that dark matter profile shapes are 
preserved in mergers. 

The most obvious comparison with observations, con- 
cerning dynamics from dark matter, is the rotation 
curves. In Fig. [¥] we plot the rotation curve for the gas 
in one of the larger galaxies (SI). The solid red line is 
the average in each radial bin of the magnitude of vec- 
tors fxv/\f\ of which there is one for each cold disk gas 
particle. The dashed lines correspond to the Keplerian 
velocities expected at each radius due to dark matter, 
baryons and the sum of the two. 

Even though the dark matter distribution in the core 
of the simulation seems to asymptote to a cusp as far 
as it is possible to ascertain above the smallest length 
scale resolved, it turns out that the baryonic mass in the 
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FIG. 4: Rotation curve of the disk gas for galaxy SI (solid 
red curve). Also plotted are the Keplerian velocities expected 
due to baryons, dark matter and in total, assuming spherical 
symmetry (v x = ^GM x (<r)/r). 



center of the galaxy is predominantly responsible for the 
rotation of the inner part of the gas disk. The rotation 
curve is comparable to observed rotation curves of large 
disk galaxies [H, [29|, [3(3] • 

The rotation curves in the simulations of the two 
smaller galaxies are somewhat too centrally peaked com- 
pared to most observations. This is probably related to 
the fact that these simulated galaxies don't fully over- 
come the angular momentum problem in the inner couple 
of kpc, and hence that the central baryonic component of 
these galaxies might still be somewhat too concentrated. 

In the inner couple of kpc there is a discrepancy be- 
tween the actual circular velocity of the gas and the to- 
tal rotational velocity expected from the enclosed mass. 
This is mainly a numerical effect, due to gravity soften- 
ing, but also partly due to effects of noncircular motions 
of, and pressure gradients in, the cold g discussed 
by P. 



IV. TESTING ADIABATIC CONTRACTION 

The most commonly used model of baryonic contrac- 
tion was suggested by Blumenthal et al. and is based 
on two assumptions, namely that the orbits of particles 
are circular and that the dark matter halo contracts adi- 
abatically, i.e. slowly, compared to the dynamical time 
scale of the system. Consider a dark matter halo which 
has an initial mass profile Mi(r). One can then ask what 
the effect would be of changing a fraction / of those par- 
ticles into baryons. The dissipational baryons will cool 
and contract, and end up with a final mass distribution 
M(,(r). In the adiabatic contraction model of Blumenthal 
et al., the relationship between the initial mass profile 



and the final dark matter profile M x (r) is given by 

r [M b (r) + M x (r)} = nM^n) = r t M x (r)/(l - /) (3) 

which relies upon the assumption of conservation of angu- 
lar momentum, a spherically symmetric gravitational po- 
tential, and the noncrossing of the circular orbits during 
the contraction, a criterion which gives (1 — /)Mi(rj) = 
M x (r). 

Since this proposal was put forward, it has been estab- 
lished that typical orbits in N-body simulations of dark 
matter halos are rather elliptical (see, e.g., [12]) so that 
M(r) changes around the orbit and M(r)r is no longer 
an adiabatic invariant. It has therefore been pointed out 
by Gnedin et al. [HJ that the relation in Eq. ([3]) could be 
modified to try to take this into account. In particular 
they argue that using the value of the mass within the 
average radius of a given orbit, r, should give better re- 
sults than that within the instantaneous radius r or the 
maximum radius at apogee r a . The average radius f for 
a particle is given by 



- 2 P r A 
r = — / — dr 

1 .L V r 



(4) 



where v r is the radial velocity, r p is the perihelion radius 
and T is the radial period. The ratio between r and f 
will change throughout the halo so Gnedin et al. param- 
eterize the average (f) (f averaged over the population 
of orbits at a given r) using the power law with two free 
parameters 



(r) = r 200 A ( — 

V r 200 



(5) 



For their simulations, they find values of A = 0.85 ±0.05, 
w = 0.8 ± 0.02 (Note that their definition of the virial 
radius is riso rather than r2oo as used in this paper, but 
the effect of this difference upon A is very small; Ai§o = 

^200 (?"200 Al80 j 1 ^""-) 

This power law assumption for the relationship be- 
tween r and (f) in Eq. (JSJ) turns out to be a good model 
of the orbital structure of all our halos, as illustrated by, 
e.g., the DM1 halo in Fig. [5] , where (f) is plotted as a 
function of r with a power law fit running through the 
data. We do not actually perform the integral in Eq. Q, 
but rather take the average of the radii over the 5 time 
snapshots and bin radially, and fit to the average within 
each bin, to determine A and w. To calculate the (A,w) 
parameters a % 2 fit is performed in 19 radial bins between 
r m j n and 0.1r2oo and using the standard deviation of the 
average of f in each radial bin. The results are presented 
in Table El 

In order to test the hypothesis of Gnedin et al. we first 
determine the best fit values of A and w directly from 
the orbital structure (or in other words the dark mat- 
ter "ellipticity" ) , as discussed above. We then perform 
baryonic contraction of the dark matter in the simula- 
tions without baryons by scanning over different A and w 
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r [kpc| 

FIG. 5: Points are (f) vs. r for the halo DM1. The best 
fit relation, corresponding to (A, w)=(0.74,0.69), shows that 
the power law assumption 0) is an excellent representation 
of the data. The red crosses represent the binned data and 
the smaller horizontal red lines represent the vertical variance 
on that data. The smaller figure shows in red ( darker shad- 
ing) the la (68%) confidence region whereas the green (lighter- 
shading) area is the 3o (99.7%) confidence region in the (A,w) 
plane. 



TABLE V: Values of A and w in Eq. by fitting (r) as a 
function of r. The ranges stated are the joint la intervals. 



DM halo 


^4min 


Abest 


Amax 


tUmin 


Wbest 




DM1 


0.71 


0.74 


0.79 


0.67 


0.69 


0.71 


DM2 


0.78 


0.83 


0.88 


0.69 


0.71 


0.73 


DM3 


0.72 


0.76 


0.81 


0.80 


0.81 


0.83 


DM4 


0.69 


0.74 


0.78 


0.79 


0.80 


0.82 



(see below), to test whether the values required to recon- 
struct the simulations with baryons correspond to those 
found directly from the orbital ellipticities. 

For each galaxy, we label the halo resulting in the dark- 
matter-only simulation by '£>M' so that the mass in- 
side radius r of that halo is denoted by M DM (r). The 
second simulation, including dark matter and baryons, 
we denote by 'S'. We thus label the mass profile of 
the baryons by M^(r), whereas the mass profile for the 
dark matter that will be predicted by the Gnedin et al. 
model will be denoted by M* m (r). The predicted 
density profile M^ m (r) is then compared to what is actu- 
ally found from our simulation including baryons, which 
we denote Mf m (r). To run the simulations with and 
without baryons with the same total mass it is neces- 
sary to have more dark matter in the simulations without 
baryons, a fact which we correct for using the parameter 
/ which represents the fraction of total mass in the form 
of baryons. The simulations use the value / = 0.15 as 



suggested by cosmology. The relationships between the 
total mass of the two simulations are thus 

M DM (r 200 ) ~ Aff m (r 200 ) + A4 s (r 200 ) (6) 
(l-/)M DM (r 200 ) * M| m (r 20 o). (7) 

The modified adiabatic contraction model is given by 

M DM «r" ))n = [M&«f/» + Mf «fy))] t } , (8) 

and the equation for the conservation of mass given by 

(l-/)M DM (r i )=M4(r / ). (9) 

It is now possible for a given A and w (which will set the 
(fi) dependency on r^, as well as the (ff) dependency on 
Tf) to use Eqs. JHJ and © to find M^Jj), using only 
M (r) and M^(r), in an attempt to reproduce M§L(r). 
That is, we can solve the above equations numerically to 
find Tf. Using the above model, we can hence determine 
predicted pinched dark matter halo profiles for any given 
set of A and w values. 

As discussed in the previous section, we use the rms 
dispersion of the number of particles in radial bins across 
the five time frames to estimate the uncertainties in our 
simulations. For each set of (A,w) values we perform 
the contraction of the dark-matter-only halo, and com- 
pare that with the halo from the simulation containing 
baryons. We restrict ourselves to the inner region of the 
dark matter halo where the effects of baryonic contrac- 
tion are most prominent. Choosing the linear bin steps 
to be lkpc (0.5kpc for the smaller galaxies) smooths out 
most of the random substructure density fluctuations and 
at the same time enables us to catch the overall shape 
of the density profiles. To be more precise, we do the 
X 2 analysis between the inner radius r m - ln and the outer 
radius 0.1r2oo (i-e. 21 kpc for the two larger galaxies and 
10.5 kpc for the smaller), and divide the range into 19 
linear bins (consequently we have Ndof = 19 — 2 = 17 
degrees of freedom (dof) in the fits). 

To illustrate which (A,w) values provide good recon- 
structions of the baryonically compressed dark-matter- 
only halo profiles, we show in Fig. [6] and [7] the x 2 of 
the la (68% confidence) region and 3cr (99.7% confi- 
dence) region in the A-w plane. In other words inside 
the black/red regions, x 2 is less than Xmin+2.3 and inside 
the gray/green regions x 2 is less than x^+11.8. For the 
canonical statistical scenario, one expects a best x 2 per 
degree of freedom value of around 1, which is also sim- 
ilar to what we find in our contour plots: Xmin/^dof = 
0.97,0.77,0.64(0.51) and 1.3 for halo 1, 2, 3(3-s) and 4, 
respectively. For any reasonable variation of the number 
of bins and radii ranges in the inner region these confi- 
dence levels stay rather stable and the best x 2 l^dof stay 
fairly constant. 

From the contour plots it follows that the fits for 
(A,w)=(l,l) ~ which corresponds to circular orbits and 
therefore the original model of Blumcnthal et al. - are 
significantly worse than the fits for the optimal values. 
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FIG. 6: Best fit parameters for reconstructing the baryonic 
compressed dark matter halo (SI) from its dark-matter-only 
halo (DM1). The red/black area is the la (68%) confidence 
region and the green/gray is the 3a (99.7%) confidence region. 
The (A,w) value, expected from the analysis of ellipticities 
as proposed by Gnedin et al. is marked by a cross, and the 
original model, by Blumenthal et al., by a circle. 



FIG. 8: Density profiles for galaxy halo SI and DM1. The 
blue/dot-dashed curve is pr 2 for the dark-matter-only sim- 
ulation, and the black/solid is pr 2 for the simulation with 
baryons. The green/dotted curve shows the reconstruction us- 
ing the model of Blumenthal et al. (i.e. (A, w) — (l,l)). The 
red/dotted curve with the error bars shows the best fit recon- 
struction, corresponding to (A,w)=(0.51,0.60). 




FIG. 7: Best fit parameters for reconstructing the i 
pressed dark matter halo (S2) from its dark-matter-only halo 
(DM2). The red/black area is the la (68%) confidence region 
and the green/ gray is the 3a (99.7%) confidence region. The 
(A,w) value, expected from the analysis of ellipticities as pro- 
posed by Gnedin et al., is marked by a cross, and the original 
model, by Blumenthal et al., by a circle. 



This can also be seen very clearly by plotting the con- 
tracted dM I dr profile using the Blumenthal et al. model 
and comparing this to the contraction observed with the 
best fit values. For example, for galaxy 1 the best fit value 
are (A,w) ~ (0.5,0.6) and the comparison are shown in 
Fig. E 

The values of (A,w) obtained directly from the relation 
between (f) and r in Eq. ([5]) are significantly different 



from what is found by the above procedure. This strongly 
suggests (not surprisingly), that there is more physics at 
work than can be described by a simple two-parameter 
model. 

From our procedure to determine the best reconstruc- 
tion of the contracted dark matter halo, it should be ob- 
vious that the confidence regions in Fig. [6] and [7] should 
not be interpreted as strict confidence regions for some 
correct values of A and w. Although the 68% and 99.7% 
confidence levels in the figures correspond to the cor- 
rect increase of x 2 , these confidence regions should rather 
be thought of simply as a separation between values in 
the A-w plane which produce good reconstruction of the 
baryonic contracted halo and those which give a worse 
reconstruction. Moreover, there may not even be any di- 
rect physical interpretation of the A and w found this 
way. 

We also investigated the adiabatic contraction of the 
two smaller galaxies 3 and 4. Again, we found that the 
relationship between r and (f) in the dark matter halos 
of these galaxies is modeled well by the power law in 
Eq. ([5]). For these simulations the best value of (A,w) 
for contraction are substantially larger, and the preferred 
A is actually even larger than those predicted from the 
ellipticities. To elaborate further on this point an extra 
simulation of the smaller galaxy S3 was performed where 
stronger early stellar energy feedback was implemented 
at a level comparable with the top of the range considered 
in [20(). This extra simulation we label S3-s. 

It was found that increasing the feedback strength (and 
therefore obtaining a less massive and somewhat more 
extended central galaxy) did change the best fit values 
of (A,w) for baryonic contraction considerably, as seen 
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FIG. 9: Best fit parameters for reconstructing the baryon 
compressed dark matter halo from its dark-matter-only sim- 
ulation, using two different types of stellar feedback. Upper 
panel is using simulation S3, whereas the lower is for an 
identical simulation, but with stronger stellar feedback (S3-s). 
The red/black area is the la (68%) confidence region and the 
green/gray is the 3a (99.7%) confidence region. The stronger 
feedback (producing a less concentrated galaxy) results in a 
significant change of the best fit values towards smaller A. 

in Fig. [HI This suggests that the details of the feedback, 
and its effect upon the concentration of the baryons, is an 
important ingredient for predicting the relationship be- 
tween the final baryonic and dark matter density profiles. 
Hence, we propose that it is not only the orbital structure 
of the dark matter halo together with the final baryonic 
profile which determines the final contracted profile of 
the dark matter halo. 

Summarizing, it is found that the parametric approach 
with different (A,w) in Eq. ©, © and © is able to very 
well reproduce the pinched dark matter profiles. How- 
ever, the (A,w) values are not universal, and do not in 
general coincide with those predicted from the ellipticity 
of the dark matter in Eq. (TJJ and ([5]). Moreover, the de- 
tails of the stellar feedback can change which values of A 
and w are preferred in the modified adiabatic compres- 
sion model. 



V. INDIRECT DARK MATTER DETECTION 

One immediate application of these results is the effect 
upon the expected indirect signal from dark matter in 
the form of weakly interacting massive particles (WIMP) 
annihilating in the galactic center [33| . Several authors 
have tried to take into account the effect of modified con- 
traction models such as the one proposed by Gnedin et 
al. upon the expected number of annihilations from the 
galactic center region (see e.g. [H, |H, [36[). It is today 
impossible for galaxy size simulations to get anywhere 
near the length resolution corresponding to the very cen- 
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FIG. 10: Diagram showing the contraction of a standard NFW 
dark matter density profile (dot-dashed blue line) by a baryon 
profile (dotted light green line) as described in the text. The 
resulting dark matter profile (dashed red lines) are plotted for 
{A,w) — (1,1) and (0.51,0.6), each splitting into two at low 
radii, the denser corresponding to the density profile achieved 
from a baryon profile which includes a central black hole. The 
extrapolated density profile from Table\T]\for the dark matter- 
in simulation SI is shown for comparison (solid black line). 
Also shown are the maximum density line and the radius cor- 
responding to the lowest stable orbit around the central black 
hole. 



ter of the galaxy. Nevertheless, we proceed in the spirit of 
comparison with the existing literature by extrapolating 
our results into extremely small radii. 

We perform baryonic contraction of two different ini- 
tial dark matter profiles with a semircalistic spiral galaxy 
baryon profile, taking some typical parameters from the 
Milky Way. To model the Milky Way baryon density 
we assume cylindrical symmetry, ignoring the possibil- 
ity of any bar. For the central bulge of stars we assume 
a density of the form p cx r~ 7 e -r / A while for the disk 
we assume a Kuzmin profile. The Kuzmin disk can be 
thought of as a delta function of matter in the z direction 
[z is the coordinate perpendicular to the disk) with a sur- 
face density OdiskM = cMdi3k °° g , where M diskoo is the 

27r(r 2 +c 2 )2 

total mass of the disk. We choose the parameters of the 
model to match observations of the Milky Way: 7 = 1.85, 
A = lkpc, c = 5kpc and with the total disk and bulge 
mass Mdiskoo = 5M bu i g0 = 6.5 x 10 10 M Q [53, 0, M, Efl] . 

The first dark matter profile we choose to contract is 
the standard NFW profile, i.e. (a,/?, 7) = (1,3,1) with 
a scale radius of 20 kpc and a local density (i.e. at r—8.0 
kpc) of 0.3 GeV cm- 3 . 

We use the baryon profile described above to contract 
the dark matter profile for some of the different values of 
(A,w) found earlier in the paper and then calculate how 
the expected flux from dark matter annihilations change. 
The results of these contractions can be seen in Fig. [TU1 
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FIG. 11: The same as Fia. \l(A but where the initial dark mat- 
ter profile instead is the exponential profile with a continuously 
decreasing logarithmic slope. The initial profile parameters 
are from the best fit to the halo simulation DM1 as given in 
TableUU 



which show both the result if a 2.6 x 10 6 M central 
supermassive black hole is included in the baryonic profile 
and if it is not. 

The same procedure is performed when the initial dark 
matter profile instead is an exponential profile, where the 
logarithmic slope is becoming continuously shallower, i.e. 
the profile in Eq. We here use the best fit values 
found for the DM1 halo as given in Table IIIII and again 
normalize the local density to be pi ~ 0.3 GeV cm~ 3 . 
The results are presented in Fig. 1111 

We do not attempt to model the complicated dynam- 
ics at subparsec scales of the galaxy [|l|, |42|, |43[ other 
than trying to take into account the maximum density 
due to self-annihilations. In other words, a galactic dark 
matter halo which has survived unperturbed by mergers 
or collisions for a time scale r ga i can not contain stable 
regions with densities larger than p max ~ mj m / (crv}T ga i- 
We assume r ga i = 5 x 10 9 years and for the WIMP prop- 
erties we adopt a dark matter mass of rridm = 1 TeV and 
an annihilation cross section of (av) = 3 x 10 _26 cm 3 s _1 . 

In Table I VII we have looked at how the different val- 
ues of (A,w) found in the previous section would affect 
the dark matter signal from the center of our idealized 
Milky Way. The Blumenthal et al. estimate gives fluxes 
far in excess of the other estimates. The other three sets 
of (A,w) in Table [VI] are, respectively, those typical of 
the best fit values found to reconstruct the larger galax- 
ies, the best fit values obtained when reconstructing the 
smaller S3-s galaxy, and the best fit to the r vs (f) curves 
which represent the ellipticities of the orbits in halo DM1. 

All our values of (A, w), when used to contract the 
dark matter halos using our analytical baryonic density 
profile, give large enhancements compared to the emis- 
sion expected from the usual NFW profile. This is true 



TABLE VI: Estimated luminosity in erg s -1 from dark mat- 
ter annihilations in the center of contracted galaxies for dif- 
ferent contraction parameters (A,w) from an initial NFW or 
the DM1 (in Table UWl dark matter profile and a baryon pro- 
file including the super massive black hole, as described in the 
text. We quota values for the inner 10 pc and 100 pc. Note 
this is the total luminosity and not that in the form of some 
specific particle species such as photons. 





NFW 


DM1 


A w 


Liopc 


Lioopc 


Liopc 


Lioopc 


Initial profile 3.9 x 10 33 


3.9 x 10 34 


2.7 x 10 31 


4.7 x 10 33 


1 1 


2.8 x 10 40 


2.8 x 10 40 


9.8 x 10 39 


9.8 x 10 39 


0.51 0.6 


2.8 x 10 37 


3.5 x 10 37 


7.9 x 10 35 


3.5 x 10 36 


0.85 0.8 


5.9 x 10 38 


6.0 x 10 38 


2.6 x 10 37 


3.3 x 10 37 


0.74 0.69 


6.9 x 10 37 


7.8 x 10 37 


1.5 x 10 36 


5.0 x 10 36 



even when the initial dark matter distribution follows 
the exponential profile as in Eq. ([2]) which does not ini- 
tially posses a cusp. From the Table [VTI we deduce that 
the boost of the luminosity compared to the standard 
NFW profile takes values in the range 10 2 to 10 4 for our 
(A, w)=(0.51,0.6) depending on the initial profile. The 
reader should however note that these extrapolations to 
very small radii neglect extra effects such as the scatter- 
ing of dark matter particles on stars or a noncentralized 
supermassive black hole [ill l42l [43| . 

The total luminosity does depend upon the (A,w) 
value, showing that the flux which can be expected 
from dark matter annihilation depends upon the extra 
physics included in simulations containing baryonic hy- 
drodynamics. 

VI. NONSPHERICITY 

In this section we relax the assumption of spherical 
symmetry and determine the triaxial structure of the 
dark matter halos, assuming ellipsoidal symmetry. 

It is straightforward to obtain the moment of inertia 
tensor 

Iij = X! i r l S ij - r i,krj,k) m k (10) 

k 

of the dark matter, gas or stars, by summing over the 
masses (m^) and positions (rk) of a matter component 
inside a spherical (or ellipsoidal) shell of a given (major) 
radius R. By diagonalizing Jy we find both the orien- 
tation of the semiaxes and the its three eigenvalues Ii. 
The three principal axes (a, 6, c) are then found from the 
relation 

a 2 - Sr ■ {-la + h + 4) (11) 

and cyclic permutations thereof, where fn is a constant 
the precise value of which depends upon the radial pro- 
file. To obtain an accurate result we then repeat the 
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FIG. 12: The triaxial parameters e = 1 — 6/a (blue/dashed 
line) and f = 1 — c/a (red/solid line) of the dark matter halos 
for the two larger galaxies with baryons (labeled SI & S2) and 
without baryons (labeled DM1 & DM2). For a perfect oblate 
shape, e = and f > whereas for a perfect prolate shape 
e = />0. 



FIG. 13: The triaxial parameters e = 1 — b/a (blue dashed 
line) and f = 1 — c/a (red solid line) of the dark matter halos 
for the two smaller galaxies with baryons (labeled S3& S4 ) and 
without baryons (labeled DM3 & DM4). For a perfect oblate 
shape, e — and f > whereas for a perfect prolate shape 
e = />0. 



search for the three axes iteratively, in each step using 
only particles inside an elliptical shell with semiaxes R, 
(b/a)R and (c/a)R as given by the previous step. Further 
iterative steps are then taken until the results converge. 
The principal axes are ordered such that a > b > c. 

Having obtained the semiaxes, one determines whether 
a halo is prolate, in other words shaped like a rugby ball 
(or an American football), or oblate, i.e. flattened like a 
Frisbee, by introducing the parameters e = 1 — b/a and 
/ = 1 — c/a. If the halo is oblate or flat it means that a 
and b are of similar size and much larger than c and the 
measure T < 0.5 where T is defined as 



Similarly, if T > 0.5 then the halo is prolate. 

The quantities e and / for the dark matter halos are 
shown in Fig. 1121 and [T3l (as a function of the major axis 
R) . It is clear that the halos change from being somewhat 
prolate (e ~ / > 0) to being somewhat oblate (e ~ 0, 
/ > 0), when baryons are included in the simulations. 
The values of T for the dark matter halos with and with- 
out baryons are listed in Table EED In the presence of 
baryons the dark matter halo tends to form an oblate 
halo, whereas when there are no baryons the dark mat- 
ter halo is more prolate. 

Given these results, the obvious thing to check is 
whether the principal axes of the dark matter distribu- 
tions and the baryon distributions are aligned. Figure [Til 
shows the alignment between the stellar disk, the gaseous 



TABLE VII: Values of the oblate/prolate-parameter T inside 
10 kpc for the halo simulations with and without baryons. 



Galaxy 


1 


2 


3 


4 


Sim. with baryons 


0.076 


0.15 


0.048 


0.24 


Sim. without baryons 


0.74 


0.92 


0.48 


0.78 



disk and the dark matter "disk" . The diagram is ob- 
tained by finding the orientation vectors of the minor axis 
(c) for different lengths scales (i.e. for different values of 
the major axis R of the ellipsoidal used in the calcula- 
tion of the moment of inertia tensor ly ) . The parameter 
A6 is then the angle between each of these vectors and 
a reference vector defined to correspond to the orienta- 
tion vector of the gaseous disk determined inside R = 10 
kpc. From the diagrams it follows that for galaxies 1, 3 
and 4, the orientation of the minor axis of the gas, stars 
and dark matter are strongly correlated with each other. 
However, for galaxy 2 the there is no clear alignment be- 
tween the different matter components at all (not even 
within the gas itself at different radii). The reason for 
this discrepancy is likely due to that the baryonic galaxy 
2 has experienced a late-time merger, making it irregular. 

Figures [T21 [JjJ] and [JJ] provide interesting information 
on the interaction between baryons and dark matter in 
spiral galaxies. Figures [12] and [TjJ] indicate that the for- 
mation of a disk by gas cooling and contraction causes the 
dark matter halo to lose most of its prolateness and in- 
stead become oblate, flattened slightly in the disk plane. 
Inside a few kpc, the dark matter can tend to be less 
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FIG. 14: Diagram showing angular alignment of the gas, stars and dark matter in our four galaxy simulations. The vertical 
scale is the difference in angle between the orientation of the smallest axis (around which the moment of inertia is the greatest) 
of the component in question relative to the axis of the gas inside 10 kpc (by definition zero and marked with a blue cross). The 
blue dotted line is the gas, the magenta dot-dashed line are the stars and the red solid line is the dark matter in the simulation 
with baryons. The black dashed line is the dark matter in the simulation without baryons, showing that the baryonic disk is 
aligned with the plane of the original density distribution. 



oblate and instead develop a prolate structure. This is 
the case for, e.g., the dark matter halo in simulation S3, 
in which it is aligned with a strong stellar bar (see [44| 
for a dedicated study of central bar structures in the cold 
dark matter). 

Additionally, Fig. [TJ] shows us that the orientation of 
the baryonic disk is rather correlated with the orientation 
of the flattest part of the dark matter halo in the simula- 
tion without baryons. The dark matter therefore seems 
to have a role in determining the orientation of the bary- 
onic disk. Subsequently, the formation and contraction of 
the baryonic disk causes most of the halo triaxiality to be 
erased, resulting in a somewhat flattened, approximately 
cylindrically symmetric halo [45l . |46( . 

The amount of triaxiality of dark matter halos is a 
generic prediction in the hierarchial, cold dark matter 
model of structure formation, and observational probes 



of halo shapes are therefore a fundamental test of this 
model. Unfortunately, observational determination of 
halo shapes is a difficult task, and only coarse con- 
straints exist. Probes of the Milky Way halo indicate 
that it should be rather spherical with / < 0.2 and 
that an oblate structure of / ~ 0.2 might be prefer- 
able (see, e.g., [47j and references therein). Milky 
Way sized halos formed in dissipationless simulations 
are usually predicted to be considerably more triax- 
ial and prolate, alth oug h a large scatter is expected 
0, HI M [Hi H M, H3. Including dissipational 
baryons into the numerical simulation, and thereby con- 
verting the halo prolateness into a slightly oblate halo, 
might turn out to be essential to produce good agreement 
with observations. In a similar numerical study [551 ] the 
effect of the baryons on the halo shape was also found to 
washout the triaxiality. However, in their one realization 
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TABLE VIII: Best fit parameters for the ellipsoidal triaxial 
dark matter halo from simulation SI (including baryons), us- 
ing model Eq. {7]). 

f e r s [kpc] a 7 xLf [46 dof] 

0.31 0.02 44.6 2.02 3.17 1.86 1.2 



of a Milky Way size halo no clear oblateness of the galac- 
tic halo was recognized, but rather an almost spherical 
halo was achieved (with / ~ e ~ 0.1). 

Having determined the ellipsoidal triaxiality of the 
dark matter distribution, one can include this informa- 
tion in the profile fits. To exemplify this, we refit the 
halo of SI (the galaxy found to be most similar to the 
Milky Way) to the 'a/37' profile as given in Eq. flU, but 
now with the replacement 

r ~* f = \j x2 + (T^P + Q?W ■ (13) 

The best fit values are given in Table IVTTT1 where we have 
used axis ratios as found within R=10 kpc (from Fig.H^l 
we note that we are quite insensitive to the exact choice 
of R). 

The (a, (3, 7) parameters of the ellipsoidal triaxial fit do 
not change much compared to the spherically symmetric 
fit, this is because the flattening of the dark matter halo 
is not very strong. However, the important difference is 
that the amount of flattening of the dark matter halo in 
the galactic plane is taken into account, and in fact, the 
goodness of the fit is slightly improved. 

The oblate structure of the dark matter will also have 
some effects on the expected indirect dark matter sig- 
nal 35]. However, the baryonic pinching effect does not 
produce such highly flattened halo profiles as the one 
proposed in, e.g., [56| to explain the EGRET observed 
diffuse gamma excess by WIMP annihilations (see also 
the critique in [57] on the dark matter interpretation of 
the EGRET signal). 



VII. SUMMARY AND CONCLUSIONS 

We have presented results comparing the structure of 
galactic dark matter halos formed in N-body simulations 
including only dark matter to that of the same halos 
formed in N-body/hydrodynamical simulations of galax- 
ies, containing dark matter, stars and gas. From our 
selected high resolution galaxies three out of the four 
galaxies formed in the hydrodynamical simulations con- 
tain very distinct disks of gas and stars, and central stel- 
lar bulges, the fourth is strongly barred. 

The central slope of the dark matter density pro- 
files becomes significantly steeper when baryons are 
present, with the average logarithmic slope parameter 



7 = — dlogp/dlogr increasing from 1.3±0.2 to 1.9±0.2 
at about 1% of the virial radius. 

The pinching of dark matter halos in response to the 
cooling and contraction of baryons was investigated fur- 
ther, to test adiabatic contraction models for the case 
where galaxies of realistic linear sizes and other proper- 
ties are formed. 

In relation to the orbital structure of dark-matter-only 
halos, the mean of the time averaged radius (f } of dark 
matter particles versus radius r is very well described 
by a power law relation (specified by two parameters A 
and w) as suggested by Gnedin et al. Moreover, it is 
found that the Gnedin et al. [12j prescription for adia- 
batic contraction is much more successful at reproducing 
the density profile of dark matter in the simulations with 
baryons, than the standard scenario of Blumenthal et al. 

in which circular orbits are assumed. However, the 
parameters of the Gnedin et al. model, (A,w), which give 
the best fit for the baryonic contraction, are somewhat 
different from the (A,w) parameters describing the av- 
eraged orbital structure of the dark-matter-only halos. 
Given the (A,w) uncertainty estimates, described in the 
text, this difference appears to be significant. In addi- 
tion, it is also found that the contraction reconstruction 
values of (A,w) also depend on the strength of the stellar 
feedback in simulations of otherwise identical halos, fur- 
ther indicating (perhaps as one might expect) that the 
effects of baryonic pinching are more complicated, than 
what can be captured in this two-parameter model. 

Our results indicate that the amount of baryonic pinch- 
ing of the dark matter halos are overestimated in earlier 
works applying the adiabatic compression model by Blu- 
menthal et al., at least for Milky Way sized disk galaxies. 
This has ramifications for predictions of the (putative) 
dark matter annihilation flux from the galactic center. It 
is found, that the flux can be reduced by several orders 
of magnitude, although baryonic contraction still boosts 
the signal significantly above the value one would expect 
on the basis of simulations containing only dark matter. 

Finally, we have determined the triaxiality of the dark 
matter halos. Dark matter only halos were found to be 
significantly more prolate than halos containing baryons. 
The influence from baryons flatten the dark matter into 
a slightly oblate halo (c/a = 0.73 ± 0.11) aligned in the 
same plane as the stellar/gaseous disk. Moreover, in 
the simulations containing baryons, galactic disks tend 
to form in the planes aligned with the flattest parts of 
the dark matter halos formed in the corresponding sim- 
ulation without baryons. 

Acknowledgments We are extremely grateful for 
conversations with Lars Bergstrom, Anne Green, Ana- 
toly Klypin, Edvard Mortsell, George Rhee and Christian 
Walck. All computations reported in this paper were 
performed on the IBM SP4 and SGI Itanium II facili- 
ties provided by Danish Center for Scientific Computing 
(DCSC). MF receives funding from the Swedish Research 
Council (Vetenskapsradet) . The Dark Cosmology Centre 
is funded by the DNRF. 



15 



J. F. Navarro, C. S. Frenk and S. D. M. White, "The 

Structure of Cold Dark Matter Halos", Astrophys. J. 

462, 563 (1996) [arXiv:astro-ph /9508025] . [19 

J. F. Navarro et al., "The Inner Structure of Lamb- 

daCDM Halos III: Universality and Asymptotic Slopes" , 

Mon. Not. Roy. Astr on. Soc. 349, 1039 (2004) [20 

[arXiv:astro-ph/031123l"1 . 

B! Moore, T~. Quinn, F. Governato, J. Stadel and 

G. Lake, "Cold collapse and the core catastro- [21 

phe", Mon. Not. Roy. A stron. Soc. 310, 1147 (1999) 

|arXiv:astro-ph/9903164] . 

J. S. Bullock et al., "Profiles of dark haloes: evolution, [22 
scatter, and environment", Mon. Not. Roy. Astron. Soc. 
321 (2001) 559 [arXiv:astro-ph/9908159| . 
T. Fukushige, A. Kawai and J. Makino, "Structure of 
Dark Matter Halos From Hierarchical Clustering. III. [23 
Shallo wing of The Inner Cusp" , Astrophys. J. 606, 625 

(2004) [arXiv:astro-ph/0306203] . [24 

D. Reed et al., "Evolution of the Density Profiles of Dark 
Matter Haloes", Mon. Not. Ro y. Astron. Soc. 357, 82 

(2005) [arXiv:astro-ph/0312544j . 

O. J. Eggen, D. Lynden-Bell and A. R. Sandage, "Evi- [25 
dence from the motions of old stars that the galaxy col- 
lapsed", Astrophys. J. 136, 748 (1962). 

P. Young "Numerical models of star clusters with a cen- [26 
tral black hole. I - Adiabatic models" Astrophys. J. 242, 
1232 (1980) 

Y. B. Zeldovich, A. A. Klypin, M. Y. Khlopov and 

V. M. Chechetkin, "Astrophysical Constraints On The [27 

Mass Of Heavy Stable Neutral Leptons. (In Russian)", 

Yad. Fiz. 31, 1286 (1980) [Soviet J. Nucl. Phys., 31, 664 

(1980)]. 

G. R. Blumenthal, S. M. Faber, R. Flores and J. R. Pri- [28 
mack, "Contraction Of Dark Matter Galactic Halos Due 
To Baryonic Infall", Astrophys. J. 301, 27 (1986). 
B.,S. Ryden, and J. E. Gunn, Astrophys. J. 318, 15 
(1987) [29 
O. Y. Gnedin, A. V. Kravtsov, A. A. Klypin and D. Na- 
gai, "Response of dark matter halos to condensation of 
baryons: cosmological simulations and improved adia- [30 
batic contraction model", Astrophys. J. 616, 16 (2004) 
|arXiv:astro-ph/0406247] . [31 
V. F. Cardone and M. Sereno, "Modelling the Milky 
Way through adiabatic compression of cold dark 
matter halo", Astron Astrophys. 438, 545 (2005) 
|arXiv:astro-ph /0501567] ■ [32 
J. A. Sellwood and S. S. McGaugh, "The Compression 
of Dark Matter Halos by Baryonic Infall" , Astrophys. J. 
634, 70 (2005) [arXiv:astro-ph/0507589] . 

E. Vasiliev, "A simple analytical model for dark mat- [33 
ter phase-space distribution and adiabatic contraction", 
JETP Letters 84, No.2 (2006) [arXiv:astro-ph /0601669] . 
A. A. Dutton, F. C. van den Bosch, A. Dekel and [34 
S. Courteau, "A Revised Model for the Formation of Disk 
Galaxies: Quiet History, Low Spin, and Dark-Halo Ex- 
pansion" , arXiv:astro-ph/0604553 

J. Barnes & S. D. M. White, "The response of a spheroid [35 
to a disc field or were bulges ever ellipticals?", Mon. Not. 
Roy. Astron. Soc. 211, 753 (1984). 

R. Jesseit, T. Naab and A. Burkert, "The Valid- [36 
ity of the Adiabatic Contraction Approximation for 



Dark Matter Halos," A strophys. J. 571, L89 (2002) 
[arXiv:astro-ph /0204164] . 

V. Belokurov et al., "The Field of Streams: Sagittar- 
ius and its Siblings", A strophys. J. 642, L137 (2006) 
[arXiv:astro-ph /0605025] . 

J. Sommer-Larsen, M. Gotz and L. Portinari, "Galaxy 
Formation: CDM, Feedback and the Hubble Sequence", 
Astrophys. J. 596, 47 (2003) [arXiv:astro-ph/020 4366 . 
J. Sommer-Larsen, "Where are the "Missing" 
Galactic Baryons?", A strophys. J. 644, LI (2006) 
[arXiv:astro-ph /0602595] . 

L. Portinari and J. Sommer-Larsen, "The Tully- 
Fisher relation and its evolution with redshift in 
cosmological simulations of disc galaxy formation", 
|arXiv: astro-ph/060653 1] 

F. Governato et al, "Forming Disk Galaxies in Lambda 
CDM Simulations" , |arXiv:astro-ph/0602351~| 

C. Power et al, "The Inner Structure of Lamb- 
daCDM Halos I: A Numerical Convergence Study", 
Mon. Not. Roy. Ast ron. Soc. 338, 14 (2003) 
[arXiv:astro-ph/0201544] . 

A. Klypin, A.V. Kravtsov, J.S. Bullok and J.R. Primack, 
"Resolving the Structure of Cold Dark Matter Halos", 
Astrophys. J. 554, 903 (2001) 

J. Diemand, M. Zemp, B. Moore, J. Stadel and M. Car- 
ollo, "Cusps in CDM halos: The density profile of a bil- 
lion particle halo", Mon. Not. Roy. Astron. Soc. 364, 665 
(2005) [arXiv:astro-ph/0504215l . 

S. Kazantzidis, A. R. Zentner and A. V. Kravtsov, 
"The Robustness of Dark Matter Density Profiles in 
Dissipationless Mergers", Astrophys. J. 641, 647 (2006) 
[arXiv:astro ph/0510583] . 

M. Persic, P. Salucci and F. Stel, "The Universal rota- 
tion curve of spiral galaxies: 1. The Dark matter con- 
nection", Mon. Not. Roy . Astron. Soc. 281, 27 (1996) 
[arXiv:astro-ph /9506004] . 

Y. Sofue and V. Rubin, "Rotation Curves of Spiral 
Galaxies", Ann. Rev. Ast ron. Astrophys. 39, 137 (2001) 
[arXiv:astro-ph /0010594] . 

E. Hayashi et al, "Disk Galaxy Rotation Curves in Tri- 
axial CDM Halos" , [ajXiv:astro-ph/0408132| 

O. Valenzuela, G. Rhee, A. Klypin, F. Governato, 

G. Stinson, T. Quinn and J. Wadsley, "Is there Evidence 
for Flat Cores in the Halos of Dwarf Galaxies?: The Case 
of NGC 3109 and NGC 6822" , |arXTv:astro-ph/0509644| 
S. Ghigna, B. Moore, F. Governato, G. Lake, T. Quinn 
and J. Stadel, "Dark Matter Halos within Clus- 
ters", Mon. Not. Roy. Astron. Soc. 300, 146 (1998) 
[arXiv:astro-ph /9801 192] . 

L. Bergstrom, "Non-baryonic dark matter: Observa- 
tional evidence and detection methods," Rept. Prog. 
Phys. 63, 793 (2000) [arXiv:hep-ph /0002 1 26| . 

F. Prada, A. Klypin, J. Flix, M. Martinez and E. Si- 
monneau, "Dark Matter Annihilation in the Milky Way 
Galaxy: Effects of Baryonic Compression", Phys. Rev. 
Lett. 93, 241301 (2004) [arXiv:astro -ph/0401512] . 

E. Athanassoula, F. S. Ling and E. Nezri, "Halo geometry 
and dark matter annihilation signal", Phys. Rev. D 72, 
083503 (2005) [arXiv:astro-ph/050463L| . 
Y. Mambrini, C. Munoz, E. Nezri and F. Prada, 
"Adiabatic compression and indirect detection of su- 



16 



persymmetric dark m atter", JCAP 0601, 010 (2006) 
|arXiv:hep-ph/0506204] . 

H. Zhao, 'A Self-Consistent Dynamical Model for the [49] 
COBE Detected Gal actic Bar", Mon. Not. R . Astron. 

Soc. 283, 149 (1996) [arXiv:astro-ph/9512064] . [50] 
W. Dehnen and J. Binney, "Mass models of the Milky 
Way," Mon. Not. Roy Astron. Soc. 294, 429 (1998) 
[arXiv:astro-ph/9612059] . 

A. Klypin, H. Zhao and R. S. Somerville, "LCDM- [51] 
based models for the Milky Way and M31 I: Dy- 
namical Models" , Ast rophys. J. 573 597 (2002). 
|arXiv:astro-ph/01 10390] 

S. M. Kent, T. M. Dame and G. Fazio, "Galactic struc- 
ture from the Spacelab infrared telescope. 2. Luminosity [52] 
models of the Milky Way" , Astrophys. J. 378, 131 (1991). 
P. Gondolo and J. Silk, "Dark matter annihilation at 
the galactic center", Ph ys. Rev. Lett. 83, 1719 (1999) [53] 
[arXiv:astro-ph/990639l] . 

P. Ullio, H. Zhao and M. Kamionkowski, "A Dark-Matter 
Spike a t the Galactic Center?" P hys. Rev. D 64, 043504 
(2001) [arXiv:astro-ph/010148l] . ' [54] 
G. Bertone and D. Merritt, "Time-dependent models for 
dark matter at the Galactic center", Phys. Rev. D 72, 
103502 (2005) [arXiv:astro-ph/0501555l . 

P. Colin, O. Valenzuela and A. Klypin, "Bars and Cold [55] 
Dark Matter Halos", A strophys. J. 644, 687 (2006) 
[arXiv:astro-ph/0506627] . 

J. Dubinski, "The Effect Of Dissipation On The 
Shap es Of Dark Halos" Astrophys. J. 431, 617 (1994) [56] 
[arXiv:astro-ph/930900l] . 

I. Berentzen and I. Shlosman, "Growing Live 
Disks Within Cosmologically Assembling Asym- 
metric Halos: Washing Out the Halo Prolateness" , [57] 
|arXiv:astro-ph/0603487l 

O. Gerhard, "Mass Distribution and Bulge Formation in 

the Milky Way Galaxy" , |arXiv:astro-ph/0608343| 

C. S. Frenk, S. D. M. White, M. Davis, G. Efstathiou, 



"The formation of dark halos in a universe dominated by 
cold dark matter", Astrophys. J. 327, 507 (1988). 
J. Dubinski and R. G. Carlberg, "The Structure of cold 
dark matter halos", Astrophys. J. 378 (1991) 496. 
M. S. Warren, P. J. Quinn, J. K. Salmon and 
W. H. Zurek, "Dark halos formed via dissipationless col- 
lapse: 1. shapes and alignment of angular momentum", 
Astrophys. J. 399, 405 (1992). 

P. A. Thomas, J. M. Colberg, H. M. P Couchman, 
G. P. Efstathiou, C. J. Frenk, A. R. Jenkins,A. H. Nel- 
son, R. M. Hutchings, J. A. Peacock, F. R. Pearce, 
S. D. M. White, Mon. Not. Roy. Astron. Soc. 296, 1061 
(1998). 

Y. P. Jing and Y. Suto, "Triaxial Modeling of Halo Den- 
sity Profiles with High-resolution N-body Simulations", 
Astrophys. J. 574, 538 (2002) [arXiv:astro-ph/0202064] . 
J. Bailin and M. Steinmetz, "Internal and Exter- 
nal Alignment of the Shapes and Angular Momenta 
of LCDM Halos", Ast rophys. J. 627, 647 (2005) 
[arXiv:astro-ph/0408163] . 

P. Bett, V. Eke, C. S. Frenk, A. Jenkins, J. Helly and 
J. Navarro, "The spin and shape of dark matter haloes in 
the Millennium simulation of a LambdaCDM universe" , 
|arXiv:astro-ph/060860"7| 

S. Kazantzidis, A. V. Kravtsov, A. R. Zentner, B. All- 
good, D. Nagai and B. Moore, "The Effect of Gas Cooling 
on the Shapes of Dark Matter Halos" , Astrophys. J. 611, 
L73 (2004) [arXiv:astro-ph /0405 1 89] . 
W. de Boer, C. Sander, V. Zhukov, A. V. Gladyshev and 
D. I. Kazakov, "EGRET excess of diffuse galactic gamma 
rays as tracer of dark matter", Astron. Astrophys. 444, 
51 (2005) [arXiv:astro-ph/0508617] . 

L. Bergstrom, J. Edsjo, M. Gustafsson and P. Salati, "Is 
the dark matter interpretation of the EGRET gamma ex- 
cess compatible with antiproton measurements?" , JCAP 
0605, 006 (2006) [arXiv:astro-ph/0602632 . 



